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ABSTRACT 


Ah modei is constructed of a repairable item inventory system where 
on-hand stocks are resupplied both by repair and procurement. The basic 
properties of the model are analyzed. A Poisson failure process 1s 
assumed, and steady state probabilities are found for operation of the 
mode? under <R,Q> and <S-1,S> inventory policies. Further, a compound 
Poisson failure process under an <S-1,S> inventory policy is also 
analyzed. Measures of supply performance, such as expected number of 
backorders per unit time, are derived for the system under the various 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


number of items in repair, X(t), at t=0 
number of items on order, Z(t), at t=0 
net inventory, Y(t), at t=0 


delay between item failure and establishment of a purchase 
order for a replacement; a random variable. 


number of items repaired and returned to stock in the interval 
(0,t]; a random variable. 


number of items ordered and received in the interval (0,t]; a 
random variable. 


probability density function of procurement lead time. 
probability density function of repair time. 

probability density function of batch size. 

Pagmeswincy dtsty tb ticn function of procurement Jead time, 
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probability distribution function of repair time. 


the number of items on hand (in stock minus backorders) at any 
time t; a random variable. 


time required to inspect an item; a random variable. 
Inventory Control Point 


inventory position, i.e., on hand + in repair + on order, at 
any time t; a random variabie. 


mean procurement lead time. 


time between placement of an order and receipt of material 
into stock; a random variable. 


number of failures entering repair during an interval of length t. 


number of failures which are lost to the system during an interval 
of length t. 


number of failures placed on order during an interval of length t. 


probability an item which has failed is repairable 





v(t) 
v4 Ct) 
Vo(t) 


W(t) 
X(t) 
Y(t) 
Z(t) 


probability a failed item is returned to the repair depot. 


probability an item which arrives at the repair depot can 
actually be repaired. 


probability an item which has failed is lost to the system. 
the amount of each order. 

probability a batch of failures is repairable. 

reorder level. 

mean repair time. 

time required to repair an item; a random variable. 


type of inventory policy whereby an order is placed for Q items 
whenever inventory position first reaches x < R. 


type of inventory policy wnereby an order is placed for S-x 
items whenever inventory position first reaches x < Ss. 


in-transit delay time, customer to inspection; a random variable. 


in-transit delay time, repair depot to stock depot; a random 
variable. 


number of failures per batch; a random variable. 


th batch; a random variable. 


number of failures in the 1 
number of batches of failures in an interval (s,stt]. 
number of batches of repairable failures in an interval (s,s+t]. 


number of batches which are iost to the system in an interval 
es.St |. 


number of items which fail in an interval (s,stt]. 
number of items in repair at any time t. 

net inventory (on hand and in repair) at any time t. 
number of items on order at any time t. 


rate of a Poisson process. 
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IT, INTRODUCTION 


Although inventory theory was one of the first areas to which 
operations research analysts turned their attention, most of the work 
has been concerned with the control of consumable items, i.e. items 
that are consumed in use. Until fairly recently the control of 
repairable items did not excite a great deal of effort. However, it 
has now been realized that, at least in military applications, repair- 
able or recoverable items are a major concern of inventory managers, 
both materially and financially. Repairables often account for a 
major portion of the total dollar investment even though the number of 
items may be small compared to the number of consumable items. Further, 
because repairable items are generaliy complex assemblies, they often 
Fequire Special stcragc, nandling, packeging, and shipment methods. 


One of the major difficulties of applying inventory theory to 


ae 


consumable items is the extreme variability of item parameters, such as 
demand, procurement lead time, sources of supply » etc. Further, the 
very fact of very large numbers of items complicates an already diffi- 
cult task. Repairable items, on the other hand, have had a tendency 
to minimize some of these probiems. The total number of items is much 
smaller and it seems somewhat easier to determine accurate parameters. 
For instance, a large number of repairable items may be classified 
electronic sub-components. Often, these can be further grouped and 
then treated as a group. The sources of supply and the sources of repair 
are generally completely known and limited in number. Repair lead times 
and procurement lead times may be often determined quite accurately with 


limited variances. Demand is generally more stable in the sense that if 





a given repairable item is a sub-component of a larger unit, it is 
generally easy to find out how many larger units are in use, and this 
at least puts an upper bound on demand. 

On the other hand, the control of inventories of repairable items 
TS Complicated by the fact of two “sources of supply", i.e., repair or 
procurement. The purpose of this thesis is to analyze a model of a 
repairable item inventory system under some common inventory policies 
assuming Poisson and Compound Poisson demand processes. The question 
of optimizing the system under the various policies is not taken up, 
Since inventory managers and other military decision-makers have not 


yet agreed as to the appropriate optimizing criteria. 





IT. A MODEL OF A REPAIRABLE ITEM INVENTORY SYSTEM 


The model used in this paper is a fairly close approximation of the 
actual repairable item inventory system in use in the Navy today. Almost 
every item stocked by Naval vessels and installations is assigned for 
technical and management control to some Inventory Control Point (ICP). 
Based on economic and technical considerations, the ICP will designate 
certain items as "repairable." These items are carried in stock, at 
various stock points or depots scattered throughout the Navy Supply 
System. Although for some items the ICP does not know how much of an 
item a certain stock depot has, repairable items are almost always 
“transaction reporting" items. This means that all issues, receipts, 
losses, etc. of that item at the Stock Depots are reported daily to the 
1cP, wnicn closely manages tne item at a centrai levei. 

Since items are often designated as repairable due to high procurement 
cost or critically low system stocks, the ICP is vitally interested in 
regaining all items which have failed. To facilitate this return of 
failed items, ships or stations requesting issue of repairable items are 
required to certify that a like number of failed items "has been or will 
be turned in." However, due to the ship's operating position, distance 
from stock points or repair depots, human procrastination, etc. there may 
be some delay in failed items arriving at the repair shop. Some items 
may be lost in transit or may never be turned in. Those that do arrive 
at the repair depot are inspected and some are found not-fit-for-repair 
and are discarded for salvage value or cannabalized for usable sub-parts. 
Those that can be repaired are eventually returned to the stock depots 


for reissue. Meanwhile the ICP, observing some items being lost to tne 





system, may decide to procure scme new items to bring the total system 
Stock back tc some desired level. The system just described is shown 


in Figure 1. 
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Figure 1. Actual Repairable Item Inventory System 


It should be noted that in actuality there usually will be several stock 
depots and probably at least two repair depots. 

However, this system can be considerably simplified with little or no 
loss of generality if it is considered from the point of view of the ICP. 
Typically, an ICP will know only when items are issued, when items are 
received into stock from repair or procurement, and when failures are 
turned in to the repair depots. From this information, the total number 
in system stock or in repair at any time can be calculated easily. The 


ICP controls the procurement of new items and thus knows the number on 
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order at any time. From these remarks, it can be seen that the various 
delays in the system, as seen from the ICP, are merely parts of either 
delay in repair or delay while on order. Therefore, it will be assumed 
that items which are repairable undergo a "repair time" which is the time 
from issue of the replacement item out of stock to receipt of the corre- 
sponding repaired item back into stock. Repair time in Figure 1 is 


mmererore 7, + [ + i + T,. Likewise, "procurement lead time" is the 


] 
time between issue and return to stock of a replacement for items which 
are lost or not repairable and thus must be purchased. 

The basic quality of failed items is whether they are repairable or 
must be replaced. Assume failed items are turned into the repair depot 
with probability Py and of those an item is repairable with probability 
Pos then failed items end up in repair with probability PyP5 and end up 


tno onoauw neacisvamant with nenahahititv f1 2. as en (1 .o.n )= ty 6. nen. 
e@ebhe Gu wase = — wh wee thea ee ow ee oe whe re ww me we ® « “os a * a ne F, nr ek ee a 4 - A. 


h | : ; Dee 
Now, to simplify the model, assume that 1): Tos I, and D are equal to 
zero, or equivalently that the ICP has complete, instantaneous knowledge 
of the location and condition of all items, and can develop repair times 
and procurement lead times as given above. 

Let mean repair time be denoted by Rae Mean procurement lead time 
will be denoted by L. Assume that Ra has a probability density function 


fp and probability distribution function F Likewise L has density 


R° 
f and distribution Fi Let PyP> = P and | - PyPo = 4- Then the systen 
described above can now be modeled as shown in Figure 2. Assume that 
failure of any item is independent of the failure or non-failure of other 
items. Further assume that the number of customer requisitions in any 


period of time is a Poisson random variable with mean rate 4. The number 


of items requested per requisition may be one each, in which case the 


TT 
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Figure 2. Model of the Repairable Item Inventory System 


process is a simple Poisson process. If more than one item can be 


~ 


areca GYOCESS 1S a COMpound PoIsgon precess. The number of 
jtems in repair at any time t is X(t). The number of items on order at 
any time t is Z(t). The number on hand, which is the number in stock 
minus the number of backorders, at any time t is denoted H(t). The 

net inventory of the system, defined as the number in stock plus in 


repair minus backorders at any time t is Y(t) = X(t) + H(t). 
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TIT. ANALYSIS OF THE MODEL UNDER POISSON DEMAND 


A. GENERAL ANALYSIS OF THE MODEL 

The number of failures in any interval (s,stt] is assumed to be 
Poisson distributed with mean At. From the assumptions in Section II, 
it is seen that every failure results in either an additional item under- 
going repair or an additional item on order. Assuming that the 
classification of failures are mutually independent and independent of 
the process, it has been shown [Parzen 1962] that the number of items 
entering repair in the interval (s.,stt] is Poisson with mean rate Ap. 
Likewise, the number of failures which will be lost to the system in the 
interval (s,stt] is Poisson with mean rate Aq. Further, the conditional 
distribution of the number of items entering repair in an interval 
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with parameters n and p. 

Now it can be shown that the number of failures entering repair in a 
given interval is independent of the number of failures which are lost 
to the system. Let N, (t) be the number of failures entering repair 
during an interval of length t and No(t) be the number of failures which 
are lost to the system during the same interval. Then the joint proba- 
bility mass function for both types of failures is: 


P(N, (4) =e oN (t) = n=x) 


P(N, (t) = x, N(t) =n) 
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xe (n-x): 


P(N, (t )=x) P(N, (t)=n-x) 


Thus Ny (t) and N,(t) are independent. 


At this point, before further analysis can be done, the processes 
of repair and procurement must be given some sort of formal definition. 
The number of items in the repair depot can be viewed as the number of 
busy servers in a queue, specifically an infinite server queue. For 
Fenda mec. the ides of ar infinite seryer queue is not 
completely satisfactory. Certainly any repair depot set up to handle a 
continuing stream of failures will probably have the ability to work on 
more than one item at a time. However, depending on the Poisson failure 
process and the repair time distribution, the probability of large 
numbers of items entering or being in repair at any time may not be 
insignificant. Wevertheless, the assumption wiil be made that the 
repair process can be considered an infinite server queue, i.e. the 
repair shop immediately starts work on all items as soon as they are 
received. 

The other specification necessary for analysis is the form of the 
reorder policy. The inventory systems considered in this paper are all 
continuous review systems, that is, all transactions are recorded imme- 


diately, and thus, the inventory manager knows the exact value of stock 
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on hand, on order, or in repair at any time. A common continuous review 
policy is the <R,Q> policy, which specifies that a constant amount Q be 
ordered whenever the inventory position reaches R. The inventory posi- 
tion 1s defined as on hand plus on order plus in repair minus backorders. 
Under this policy the procurement process may be thought of as a queue 
where the number of busy servers at a time t represents the number of 
orders of size Q that are outstanding at that time. Another frequently 
used policy is the <s,S> policy, meaning order up to S whenever the 
inventory position falls below s. In the case where s = S-], an order 
is placed everytime a failure cannot be repaired. Under the <s,S> policy, 
Hiempmocunemenl Drocess 1S also comparable to an infinite server queue 
where the number of busy servers at a time t is the number of orders 


outstanding at that time. When items are demanded one at a time, the 


mm ecltcy ts the same as an <s,.S> policv. 


Now let X(t) be the number of items undergoing repair at time t and 
Z(t) be the number of items on order at time t. Then 


X(t) = N, (t) + C, - D, (t) 


Z(t) Ny (t) + C, - Dp(t) 


where C, pom Cliienc. tS) 7 (0). D, (t) is the number of items repaired and 


2 
returned to stock in the interval (0,t]}, D,(t) is the number of items 
ordered and received in the interval (0,t], and No (t) is the number of 


items lost to the system which have been placed on order by the time t. 
t 1 


N, (t) is seen to be related to N,(t) Bye 


t No(t) - 
No (t) = Q Q 


where [In] is the largest integer less than n. If items are reordered 





every time one is lost to the system, i.e. Q=1, then N, (t) = No(t). 





Otherwise No (t) < Nott). Now both X(t) and Z(t) can be likened to the 
number of busy servers in a queue. D, (t) depends only on Ni (t), the 
repair queue discipline, and repair time distribution. No (t) and D,(t) 
depend only on No(t), the order queue discipline and procurement lead 
time distribution. Assuming independence between the two queue opera- 
tions, it is seen that X(t) and Z(t) are independent. 

Let H(t) be the number of items in stock ready for issue at time t, 
where H(t) is negative if backorders exist. The net inventory, Y(t) = 
X(t) + H(t), is the number of items in stock plus in repair minus back- 
orders. It is important for further analysis, that the independence of 
X(t) and Y(t) be shown, since H(t) will be obtained from X(t) and Y(t). 
This independence has been shown by Richards [Ref. 5] and it is instruc- 
tive to repeat the proof here: 

',..8Mamine the manncr in which transitions are made in the two 
stochastic processes, {X(t), t > O} and {Y(t), t > O}. 

Whenever a type-1 failure [into repair] occurs, H(t) decreases by 
one, but X(t) also increases by one for a net change of zero in Y(t). 

On the other hand, H(t) decreases by one while X(t) remains constant when 
a type-2 failure [procurement] occurs; thus a type-2 failure causes a net 
decrease of one in Y(t) but does not affect X(t). 

The arrival of a procurement increases Y(t) by the amount Q, while 
the repair level is unaffected. Conversely, X(t) decreases while H(t) 
increases by the like amount whenever an item leaves repair. Thus, 
repairs decrease X(t) but leave Y(t) unchanged." 

Then it can be seen that 


X(t) = C) + Ny (t) - D,(t) 


Y(t) = Cz - No(t) + Do(t) 





where C, = ¥(0). As previously shown D, (t) depends only on N,(t) and 

not on N,(t). Conversely, No (t) and D,(t) are functions of No(t) but 

not N,(t). Thus X(t) is a function only of N, (t) while Y(t) is a function 
only of No(t). Since N, (t) and N,(t) have been shown to be independent, 


then X(t) and Y(t) are independent. 


B. AN <R,Q> POLICY WITH CONSTANT PROCUREMENT LEAD TIME 

Consider first the operation of the model under an <R,Q> policy, 
where the inventory position I(t) = X(t) + H(t) + Z(t) = Y(t) + Z(t). 
An order is placed for Q units whenever I(t) falls to R. Assume that 
all delay times other than repair time and procurement lead time are 
zero. Further assume that procurement lead time is a constant, L. Then 
the "net demand" placed on the system stock, Y(t), i.e. the number of 
items which must be replaced through procurement is a Poisson process of 
rate gq. Thus the model from this point of view is Merely an inventory 
System operating under an <R,Q> policy experiencing Poisson demands at 
rate gq. This model has been solved by Hadley and Whitin [Ref. 2]. 
Translating their notation into that used in this paper, the following 
are given: 


PLY(t)=y] = ~ [F(Rt1-y3AqL) - F(R+Q+1-y;AqL)] ~~ <y < Rt] 


£D|— 


[1 - F(R+Q+1-y3AqL) ] Rt] < y < R+Q 


| 


oe i-r 
where F(asr) =1- 2% {r) e 


i.e, the complementary cumulative of a Poisson distribution of rate r. 


Now 


ey = X(t) + H(t) 





thus 
Mee = Y(t) =2X(t) 


PLH(t)=h] = PLY(t)-X(t)=h] 


5 PLY(t)=ith,X(t)=i] 
i=0 


Since it has been shown that Y(t) and X(t) are independent 
PCH(t)=h] = ©  PLY(t)=ith] PEX(t)=i] 
i=0 


Note than an order is placed for Q units whenever the inventory position 
falls to R. Thus I(t) varies between R+ 1 and R+Q. But 

Peels W(t) < R + 0 
implies 


pee} < A(t) + eS R + Q 


Z(t), the number of items on order, is always non-negative, i.e. Z(t) > 0 


Thus 
tees y(t) < Rk + Q 
implies 
Mee het Q = Zit) <R +O 
or Max Y(t) = R + Q when Z(t) = 0 


Therefore P[Y(t)=y] = 0 for y> R+Q 


R+Q-! 
means PLH(t)=h} = = PLY(t)=ith] P[X(t)=i] 


1=0 
The probability that Y(t) = y has already been developed. Now P[X(t)=x] 


is merely an application of Palm's queuing theorem which states that ney 





an infinite server queue, with Poisson arrivals of rate 4 and service 
times distributed F(t) with mean M, the steady state probabilities of x 


busy servers are Poisson with rate AM, i.e. in this case, 


(xpR, )%e™*PRA 
P[X(t)=x] = -——---—— 
x 


Since mean repair time is taken to be Ras and of course 


E[X(t)] = ApR 


A 
Then 
R+Q-h 
PIH(t)=h] = x PLY(t)=ith] P[X(t)=i] 
1=0 
R+Q-h _ 
= q LF (RFT -i-hsaqh)-F(RFGHT-i-h3 ql) J PC isApRq) for ith =< Rr] 
i=0 
R+Q-h _ 
= q Ll ~ F(RtQ+1-i-h3AqL) ] p(isApRy) Ole RES eis Ran’, 
i=0 
a-r 
where p(as3r) = Keer 
a. 
Simplifying 
R+Q-h 
PLH(t)=h] = q ® [1 - F(R-h-i3;Agk) - 1 + F(R+Q-h-i3Aql)] - p(isAPRq) 
i=0 
ce , 
= q y® [F(R+Q-h-i3AqL) p(isApRy ) - F(R-n-i;AqL) p(isApRy)J 
i=0 
R+Q-h R+h . 
=a = plisdpRy) F(R+Q-h-isAql) - 2 p(isApRa) F(R-h-isAql) 
Q {=0 =O) 


Since H(t) = Y(t) - X(t) it is easy to find the expected value of stock 


on hand: 


Wy 





ECH(t)] = ELY(t)] - ELX(t)] 


Richards showed that the expected value of net inventory for an <R,Q> 
foley 1S R + aa m where m is the mean number of failures lost to 


the system in a procurement lead time. Thus 


p+ Sah - apr 


E[H(t)] ; 


R + ot APR g + gl) 


Now inventory position fluctuates between R+1 and R+Q, j.e., 
Rie Ze ete REO 


Assume the system starts with Y(0) = R + Q and thus Z(0) = 0. Then when 
Y(t) drops to R, Z(t) becomes Q. If this order does not come in by the 
time Y(t) = R - Q, then Z(t) beccmes 20 and so forth. the tablic below 


Shows this relationship: 


Ges Z(t) 
R+Q 0 
R +] 
R 
ep ey | 0 
=f 20 
R - 20 + ] 20 
R = (n-1)Q nQ 
R-(nQ) + 1 nQ 


20 





Then Z(t) = nQ if and only if R-nQtl < Y(t) < R-(n-1)Q Nelsen. 4% 


P[Z(t)=nQ] = P[R-nQt1 < Y(t) < R-nQ+Q] 


Fy (¢y (R-n@Q) Fc) (R-nQ) 
her F a op PIY(t)= 


Note that if Q=1, i.e. reorder every time an item is lost to the system, 


this reduces to 


P[Z(t)=n] = Fey (Rontl) - Fy) (R-n) 


Hl 


PLY(t)=R-nt+1] 


Hy 


g EF(nsaql) - F(nt15aqi)] 


That is, when Q=1, the number on order at any time t, like the number in 
repair, 1S Poisson aistributed as expected by Palm's Theorem. 
Now Gey = Y(t) + Z(t) 

eee ee Y(t) 1 + £[ Z(t) ] 


p+ Pe pe MO gh + Efz(t)] 


or FE[Z(t)J= AaL 


Expressions have now been developed for X(t), Y(t), Z(t) and H(t) and 
the system is thus completely described. One measure of effectiveness 
that is almost always calculated for inventory models is the expected 


number of backorders at an arbitrary point in time. It is, of course 


2] 





. b p[H(t)=-b] 
b= 


where b = the number of backorders. Under previous definition the 
quantity on hand is considered a negative value when backorders exist, 
thus H(t) = -b denotes b backorders on the records. 

Another well-known measure of effectiveness is the expected number 
of backorders per unit time which is also given by the expected number 
of failures minus the expected number of fills. The expected number 
of failures is just A, the rate of the Poisson failure process. Now 
a fill can occur only if H(t), stock on hand, is greater than zero. 

ime jes Os Y(t) - X(t) > 0 
=> V(t) > X(t) 


Now Y(t) is never greater than R+Q. If Y(t) > X(t), but X(t) > 0 
(always), then ¥(t}) > YU. So 
Dia ie—= ply(t) > X(t) for O< Y(t) < R+Q 
Therefore: 
R+Q 


+ 
pena ta) = a pLY(t) > X(t)[¥(t)=y] pLY(t)=y] 
y= 


pLX(t) < y] 
y-1 (ApR,)%e™*PRA 
= > eee a 


pLY(t) > X(t)[¥(t)=y] 


x=0 x 


iI 


F(y-13ApRq) 


Then 
R+Q 
p[fill] = : F(y-13ApR,) pLY(t)=y] 
ye 
R R+Q 
= © ply(t)=y] Fly-15apR,) + £  ply(t)=y] Fly-1sapRy) 
y= y=Rt] 
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= 2 CF(R#+Q-ysAql) - F(R-ysAql)] F(y-13ApRy) 


ye 
R+Q 
+ fF q F(R+Q-y5AqL) F(y-15ApRy ) 
y=Rt+] 
R+Q R 
=~ 2X F(RtHQ-y3;Aql) F(y-13ApR,) - © F(R-y;Aqlh) F(y-13ApR,) 
Q y=] A y= A 


Then the expected number of backorders is A- A: p[fill] 


R+Q R 
2 F(RHQ-y3AGL) F(y-13ApR,) - 2 F(R-y3Aqgl) F(y-13ApRy) 
Me y=] 


O|> 


C. AN <S-1,S> POLICY WITH STOCHASTIC REPAIR AND PROCUREMENT LEAD TIMES 
The <S-1,5> policy has often been suggested as especially appropriate 


for repairable items since the object of the policy is to maintain the 


fis. he YT S12 re ndtn ni, a i 4 i < . nm e s 
ey oe Sex a. thee level S_ 1 Ne. When an item fails and cannot. be 


ts 
oa 


Se we we 


repaired another is immediately placed on order. Assume that there is 
no lag between a non-repairable failure and the placing of an order. 


Then for this model 
yer Hit) + Z(t) = S for all t > 0 


The number in repair at any time t, X(t), and the number on order at any 
time t, Z(t) are modeled as infinite server queues, with arrival rates 
Ap and Aq respectively, where q = l-p. In this model both repair time 
and procurement lead times are random variables with means Ra and L 
respectively. Thus. Palm's queuing theorem is applicable to both X(t) 
and Z(t) giving 
(ApR, )%e7*PRA 
PEX(t)=x] = ——--—-——- Me le2 . 5.0 


X. 
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z_-AQqL 
PIZ(t)=z] = sAgbye z=0,1,2,..- 


va 


j.e. both X(t) and Z(t) are Poisson distributed. 


Now note that since X(t) + H(t) + Z(t) = S 


MiGeet Z(t) = Ss 


then = p[Y(t)=y] = plZ(t)=S-y] y an integer, -»<y<S 


7 (xq 2% AM 


(S-y): 


Also pLH(t)=h] = pLY(t)-X(t)=h] h an integer, -~<y< S&S 


S 
Xe plY(t)=y,X(t)=y-h] 
y=h 


z pLY(t)=¥1 pLXit )=y-hJ 
v= 


S Coben an we (ApR,)YNe*PRA 
», 


pLH(t)=h] [Ian — as) lam 


y=h 
Feeney and Sherbrooke [Ref. 4], using a slightly different model, 
with an <S-1,S> policy, suggested the following three measures of supply 
performance (redefined for the present model): 
(1) R(S), called the ready rate, is the steady state probability 
that S or less units are in resupply. For the present model "in resupply" 
means in repair plus on order. 
BPGt) Z(G) eS] 
S 


® pLX(t)+Z(t)=w] since -both X(t) and Z(t) are non-negative. 
w=0 


R(S) 


Si 
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I 
Bag Ue 8 2 


; p(X(t)=x,Z(t)=w-x) 


=O | x=0 
S W 
a Z plX(t)=x] p[Z(t)=w-x] 
w=0 x=0 
S ow (ApR, )*e7*PRA Gaie ees 
R(S) =x Co sean 
w=0 =x=0 ae (w-x)! 


(2) F(S), called fills, is the expected number of demands per time 
period that can be filled immediately from stock on hand. As in the 
last section this is just the failure rate times the probability of a 


mill, te. X plfill] 


pL fill] = p[H(t) > 0] 


pLX(t) < ¥(t)] 


2 pEX(t) < ylY¥(t)=y] ply (t)=yJ 
y= 


S 
Z p[X(t) < y] pLY(t)=y] 
ye 


S S-y -AgL 
r F(y-15ApRy) (agl)" “e 


Vet (S-y): 


S S-y_-AgL 
aS: F(y-15ApR,) Wel) eta 


y=] (S-y)! 


F(S) 


Now, it is well known that for a Poisson demand distribution the <R,Q> 
and <s,S> policies are the same, with R + Q = S and Q = S-s. In this 
Gaseemuien diees-1,S> policy R+Q=S, R=S-+-1, andQ= 1. In the 


previous section it was shown that for the <R,Q> model, 


Zo 





R+Q R 
z  F(R+Q-y3AqL) F(y-13ApRy) - £ F(R-y3;Aql) F(y-15ApR 
= 


pLFIN] = 5 : 
y y= 


q) 


Substituting the values for R+Q, R, and Q: 


S S-1 
p[fill] = 2 F(S-y3AqL) F(y-15ApRy ) - Y F(S-y-13AqL) F(y-13ApRaq ) 

y= y= 
S-1 

- . [F(S-y) - F(S-y-1)5Aql] F(y-15ApRy) + F(O3AgL) F(y-13ApRq) 
y= 
S-] 

= & p(S-ysAql) F(y-13ApRy) + F(O3AgL) F(y-15ApRy) 
y=] 


Now F(QO;AqL) = p(Q;AqL) = p(S-S3AqL). 
Thus 
S 
pPLFITI] = & p(S-ysAql) F(y-13ApR,) 
y=l 
which agrees with the expression given above for the <S-1,S> policy. 


The expected number of backorders in a unit time is, again, 


A- A: pLfill] 
S 
= A -d EZ Fly-135ApRy) p(S-y3Aql) 
a 
(3) A third measure of supply performance considered is S(S), the 
number of "units in service," defined as the expected number of units 
in routine resupply at a random point in time. Routine resupply is 
further defined as the number of units in resupply minus the number of 
units backordered. Again “in resupply" is taken to mean in repair plus 
on order. We note that if X(t) + Z(t) = r < S, then H(t) > 0, i.e., 


there are no backorders, and the number of units in routine resupply is r. 
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On the other hand if X(t) + Z(t) = r = S, H(t) < 0, i.e., there are 


backorders, and the number of units in routine resupply is S. 


Thus 
S 0 
S(S) = = r PLEX(t)+Z(t)=r] + © S PEX(t)+Z(t)=r] 
r=0 r=§$ 
S r ) 00 r 
= Y r2?Z PLX(t)=w] p[Z(t)=r-w]f'+ 2 S24 2m PLX(t)=w] p[Z(t)=r-w] 
r=0 w=0 r=S w=0 
S r (ApRy)e™APRA (ag)? Me™ “A> 
r=0 w=0 Ww! (r-w)! 
a r (Ap, )We7APSA (nah )o Me AML 
+ S » . o— 


r=S }w=0 Ww! (r-w)! 
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IV. ANALYSIS OF THE MODEL UNDER COMPOUND POISSON DEMAND 


The model in the last section assumed that failures were a Poisson 
process of rate 4, i.e. the number of failures in a time interval of 
length T has a Poisson probability mass function with parameter rT. 

This implied failures occurred one at a time, and thus each requisition 
for replacement was for exactly one item. 

In practice, the Poisson mocel may be quite reasonable for major 
repairable units, i.e. items of high cost and/or large, complex structure. 
However, there are smaller not so expensive items that are also economi - 
cally repairable and that may experience failures in batches, i.e. more 
than one at a time. One good example might be the printed circuit “cards” 
found in many modern computers, radars, and other electronic equipment. 
jnese cards aré generality easily removed and insta 
tions of peak load, current surges, etc., several may fail at the same 
time. Generally, the only damage may be to the resistors, capacitors, 
transistors, etc. which are soldered onto the printed circuit, while the 
circuit card itself is undamaged and repair is feasible. 

These types of items lead to the consideration of failure processes 
in which the number of failures in a given time period is a compound 


Poisson process {W(t), t>0}. This process can also be represented by 


V(t) 
W(t) = x  U 


n=] i 


where {V(t), t>0} is a Poisson process which is the number of customer 
requisitions received and {U_; Ne ieceeee iS deka ly Of independent 
identically distributed random variables where U, is the number of items 


th 


requested or the n° requisition. It is necessary that fV(t), t>0} and 
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U, are assumed to be independent and for this model that means that 
the number of items per requisition is independent of which requisition 
it is, which 1s an entirely reasonable assumption. It is noted that the 
Simple Poisson process is just a special case of a compound Poisson 
process, where the probability {U=1,for any n} = 1. The compound 
Poisson is an extremely rich family of distributions. Various seemingly 
quite different processes can be described merely by the proper choice 
Of a compounding distribution. 

Assume that W(t) the number of items which fail in a unit time 
interval is a compound Poisson process with rate A andacnee U, is 
distributed as a random variable U, with probability mass function iy 
and characteristic function dy: Then it has been shown [Parzen, p. 130] 


mat the Characteristic function of W(t) is 


i = CHC Stee 
P(E) Xe? e : 


Referring to the repairable item inventory model of this paper, 
assume that all failure items in a particular batch are either repairable 
or all are non-repairable and result in a procurement. Further, let the 
probability that all items in a batch are repairable be r. Then the 


number of items to enter repair from the nth 


batch, Ry? 1s 
U, with probability r 
Q with probability l-r 

That is prob [R, = Ol aos EU = Mul US len since 


prob [R, = 0] = I-r 
The characteristic function of R, 1s 


a= EfelMny = * el) p[R as] 
j=0 





R =O} + x _idj 
pl n ] ne ra ?J pLR =3] 
ro eae ‘il g'PJ plu =3J 

j= 


l-r + ro 
U 
n 


PR =l-rt roy 


Then N,(t), the number of items entering repair in a time interval (0,t], 


is 
(t) 
yeh 
n=] 


V 


Ny (t) n 


and from above 


fae N 
(= 
/ 


are wate a \ 

? 

Nw ry £5 A vt 
"UN ~ 


mney = ls 


2 ert(l-rtrdy-1) 


- ,dtr(dy-1) 


which 1s just the same compound Poisson process as failures with rate A 
replaced by rate Ar. Likewise, No(t)s the number of items which are 
placed on order in an interval (0,t], can be shown to be the same 
compound Poisson process with rate A(I-r). 

Now, let the number of batches entering repair be V, (t) and the 
number of batches which must be replaced by procurements be Vo(t), SO 
that V(t) = V(t) + Vo(t). The number of batches in the compound 
Poisson model corresponds to the number of items in the simple Poisson 


model. In Section II it was shown that the number of items entering 


repair, No (t), and the number of items lost to the system, No(t): were 





independent random variables. In the compound Poisson process, this 
implies that V(t) and Vp (t) are independent. Let Ny (t) be the number 
of failures in an interval (s, s + t) which can be repaired, and tet 
No(t) be the number of failures in the same interval which are los ato 


the system. Then 


eres” 
N(t) = 2 Ux 
n=1 n't 

V(t) 

No(t) = U Xo 

n=| 

where -)! if the nt” patch goes to repair 
X} “)0 otherwise 


Oo ii the ae batch goes to repair 
Xo “\1_ otherwise 


m oo Sacred Spe { 
Dy assuiietion wb 


eed 


n? n=t,Z,...} is a family of ineepenscns tdencicadly 
distributed random variables, and a batch is either repairable or lost 
to the system. Thus Ny (t) and No(t) are sums of independent identically 
distributed random variables, and, furthermore, no Uy, which is in the 
SV eels N,(t) is also in the sum for No(t) and vice versa. Therefore, 
N(t) and No(t) are independent. Further, X(t), the number of items in 
repair at any time t is a function only of N,(t) and the repair queue 
parameters. Likewise, Z(t), the number of items on order at any time t 
is a function only of No(t); the procurement policy, and the procurement 
queue parameters. Thus VGalvand 2(t) anemia independent for the 
compound Poisson case, assuming a batch goes completely to repair or is 
completely lost to the system. 


Again defining H(t) to be the number of items on hand at time t, and 


calling Y(t) = X(t) + H(t) the net inventory at time t, the same analysis 
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as in Section II shows that X(t) and Y(t) are independent. Using the 


notation of that section 


X(t) = C, + Nj(t) - D,() 


Y(t) = Cy - No (t) + D,(t) 


where C, = X(0), C. = Y(0), and Ny (t) is the number of items which have 
been ordered in the interval (0,t]. It has been shown that Nj (t) and No(t) 
are independent and it is assumed that D, (t) 1s a function only of Ni (t), 
1.€., that the number of items repaired in an interval of length t is a 
function only of the number of items which enter the repair queue and the 
queue operation itself, and that No (t) and D,(t) are functions of N,(t). 
Thus X(t) and Y(t) are functions of independent random variables and are 
themselves independent random variables, 


Re Ce tet ak Js ey me) cma os oe 
MoS lime LIS ffl re ey 


4q e 
v ! 
Ld 


* pipe Gel OCS nniarw 
=e on hd “Aw © © wy 


Ko 
~ ¥ Abe 


whereby every failure that cannot be repaired is immediately placed on 
order. This, of course, may result in an order for more than one item 
under the assumption that an entire batch is either repairable or not. 
Assume further that all items which enter repair or which are placed on 
order at the same time have the same repair time or procurement lead 

time respectively. The effect of these assumptions is that items which 
fail in the same batch are either repairable and are repaired as a batch, 
or that those which are not repairable must be ordered and are received 
into stock as a batch. While these assumptions may seem unduly restric- 
tive, in actual practice they are quite reasonable. Certainly items 


ordered at the same time are generally received at the same time. Further, 


repair often consists of a series of operations to be performed. It is 





not unusual that the first operation is completed on all items in a 
batch before the second operation is started, etc., and upon completion 
of repair the items are returned to stock as a batch. 

Let Ry be mean repair time and L be the mean procurement lead time. 
Let F Ax) be the "compounding" probability mass function, that is 
F Ax) = P{U = x}. 

Now Feeney and Sherbrooke have shown that under the above assumptions 
given a compound Poisson arrival density of rate A into a queue, the 
Steady state probability of x units in the queue is the same compound 
Poisson density with rate AT, where T is the mean service time in the 


queue. Since the rate of items entering repair is Ar: 


Prob {x units in repair at any time} 


Prob{X(t) = x} 
(avr) Te APRA 


Fn f(x) 
i=0 i ° 


1* 
where fi is the i-fold convolution of fie 


Likewise, with the rate "into" procurement being A(1-r), 


Prob{z units on order at any time} 


Prob{Z(t) = z} 


1* 


arate (sel ; 
0 i “ 


(z) 


lM N 


a] 


A 


Then assuming the <S-1,S> policy, where an item is ordered every time it 
is lost to the system, the inventory position, which is equal to on hand 
+ in repair + on order - backorders, is always equal to S. In terms of 
the model 

H(t) + X(t) + Z(t) = S 


33 





Oper a) altel == S 


Then Prob {Y(t)=y} 


Prob{S-Z(t)=ys 
= Prob{Z(t)=S-y} 


i -A(Ier)L ogg 
sy BO-r)L1e : Fi (sy) 


=O cle 


Furthermore, since Y(t) = H(t) + X(t) 
Rice Vit eee ( t ) Soret) = Ss 


Then Prob {H(t)=h} 


ProblY(t)-X(t)=h} 


S-h 
XE Probly(t)=yth,X(t)=y} 


= 


S-h 
x Prob{Y(t)=yth} Prob{X(t)=y} 
y=0 


and since Y(t) and X(t) are independent 


Sch }S-y-h [A(ler}LJie UPL oe y (Arka) te7 TR an 
= Ef XS #s-y-n) 3 f(y) 
ye | dat ie eo 


The system is now described for any time t > 0. Again, the measures of 
supply performance suggested by Feeney and Sherbrooke are ready rate, 


fills and units in service. 


Ready rate, R(S) is the steady state probability that S or less units 


Paeceienenair Of om order, i.e. X(t) + Z(t) < S. 


Prob [X(t) + Z(t) < S] 


N 


R(S) 


S 
x Prob [X(t) + Z(t) = w] 
w=0 





5 W 
2B 


= Prob [X(t)=x,Z(t)=w-x ] 

w=0 x=0 

S W 

. y Prob [X(t)=x] Prob [Z(t)=w-x] 
w=0 x=0 

SW X ute ie ix WX EtieAn lela - ‘| 
— 6S % (x) Yo ——————-———-—— fF  (w-x) 
w=0 x=0 | i=0 i. i=0 ie : 


Fills, F(S), is the expected number of demands per time period for an 
item that can be filled immediately from stock on hand. Feeney and 
Sherbrooke reasoned that, "There will always be one fill if S-1 or less 
units are in resupply [in repair plus on order] since each customer 
places an order for at least one unit. There will be a second fill if 
S-2 or less units are in resupply and tne customer order is for two or 
Pees. 

Therefore, F(S) is the expected number of customers per time period, 
i.e. 4, times the expected number of filis per customer. As shown above, 


Piiiascan) DE wepresented as 


F(S) = A f (probability a customer's demand is greater than 0) - 
(probability X(t) + Z(t) < S-1) 
+ (probability a customer's demand is greater than 1) - 


(probability X(t) + Z(t) < S-2) 


#e gee 


+ (probability a customer's damand is greater than S-1) ° 
(probability X(t) + Z(t) < 0.) 
= AST © R(S-1) + [1-F,(1)IR(S-2) + [1-F,(U)> Fy (2)]R(S-3) + 
+ [1-f (1)-#, (2)-...-f (S-1)] R(0) f 





This is the same formula given by Feeney and Sherbrooke, except the ready 
rate has been adjusted for the double “sources" of supply, i.e. repair 
and procurement. Then the expected number of backorders in a unit time 
is, with probability one, the expected number of demands less the expected 


number of fills, i.e. 
E[number of backorders] = AE[U] - F(S) 


The other measure of supply performance is units in service, S(S), which 
is the expected number of units in routine resupply at a random point in 
time. There are r items in routine resupply if X(t) + Z(t) =r < S, and 


Emene dre S items in routine resupply if X(t) + Z(t) > S. Thus 


S e 
oye > Y Probix(t)+Z(t)=r] + & S Probl X(t)+Z(t)=r] 
r=] r=S+] 
S r 7 
=- Z£ r £  ProblLx(tj=wy Probl Z(Li=r-w i 
r=] = w=0 
00 , 
+ S x Prob[X(t)=w] Prob[Z(t)=r-w] 
skal w=Q 
S ie W (APR, \\¢ AEN ax r-w [AC] -r)Ly'e “AC =r DL a 
=> + a A +! (Ww) 5 —————_F, (r-r 


0 it i=0 i 


ed Yr W (ary) 'e7 APRA {* 
+S & 4 z i aa Ui (w) ° 


i=0 1! 


36 





V. CONCLUSIONS 


The results developed in this paper are relatively straightforward 
extensions of models and results found by Hadley and Whitten and by 
Feeney and Sherbrooke. The model is a fairly good approximation of 
reality. The assumption of Poisson or compound Poisson demand processes 
has found much favor in inventory theory. The steady state probability 
of the number of items in repair, in stock, or on order, have been 
found, under various inventory policies. From these probabilities, 
formulas for several measures of supply performance have been derived. 
Other measures can be developed in a similar manner. Thus, given the 
steady state system probabilities and the derived measures Of effective- 
ness, inventory managers can now attempt to optimize the system under 
way ous criteria. 

Further analysis to improve the model may also be possible; for 
instance, the probability of a fill given a delay time ¢ in which to 
make the fill would be a useful result. It may also be possible to 
relax the assumption that under compound Poisson demand, failure batches 
either go entirely to repair, or are lost to the system. It would seem 
more natural that each item of each batch may go into repair or may be 


lost to the system. 


37 





BIBLIOGRAPHY 


Galliher, H. P., Morse, P. M., and Simond, M., "Dynamics of Two 
Classes of Continuous Review Inventory Systems," Journal of the 


Operations Research Society of America, v. /, Pp. 362-384, June 1959. 


Hadley, G. and Whitten, T. M., Analysis of Inventory Systems , 
Prentice-Hall, 1963. 


Parzen, E., Stochastic Processes, Holden-Day, 1962. 


re 


Rand Corporation Memorandum RM-4176-1-PR, the (S-1,S) Inventory 
Policy Under Compound Poisson Demand: A Theory of Recoverable Item 
Stockage, by G. J. Feeney and C. C. Sherbrooke, March 1966. 


Richards, F. R., Ph.D. Dissertation, Clemson University, SVE 


Takacs, L., Introduction to the Theory of Queues, pp. 159-162, Oxford 
University, 1962. 





INITIAL DISTRIBUTION LIST 


No. 


Chief of Naval Personnel 
Pers-116 

Department of the Navy 
Washington, D. C. 20370 


Naval Postgraduate School 
Department of Operations Research 
and Administrative Sciences 
Monterey, California 93940 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Assistant Professor F. R. Richards, Code 55 
Department of Operations Analysis 

and Administrative Sciences 

Naval rostgraduate Scnco} 


Monterey, California 93940 


LCDR Michael G. Wilson, SC, USN 
Defense Electronics Supply Center 
DESC-DD 

Dayton, Ohio 45401 


Professor P. L. Zehna, Code 55 
Department of Operations Analysis 
and Administrative Sciences 

Naval Postgraduate School 
Monterey, California 93940 


Commanding Officer 
Fleet Material Support Office, Code 97 
Mechanicsburg, Pennsylvania 


Commander 

Naval Supply Systems Command 
SUP 0613 

Washington, D. C. 


oF 


Copies 


Be 





Security Classification 
DOCUMENT CONTROL DATA-R&D 


(Security classtlication of titie, body of abstract and indexing annotation must be entered when the overall report is classified) 


1 ORIGINATING ACTIVITY (Corporate author) 2a, REPORT SECURITY CLASSIFICATION 


Unclassified 


Naval Postgraduate School 
Monterey, California 93940 


3 REPORT TITLE 


A Stochastic Repairable Item Inventory Model 


» 


DESCRIPTIVE NOTES (Type of report and inclusive dates) 


Master's Thesis; September 19/1 


§. AUTHOR(S) (First name, middie initial, last name) 


Michael George Wilson 


fe REPORT DATE —— = =. Ja. TOTAL NO. OF PAGES ) 
tonnoepeember 197 | 4] 6 : 


Ba. CONTRACT OR GRANT NO 19a. ORIGINATOR'S REPORT NUMBER(S) 





- PROJECT NO. 


o 
i Sb ee 


ms lob. OTHER REPORT NO(S) (Any other numbers that may be assigned 
this report) 


a 


} 10. DISTRIBUTION STATEMENT 
Approved for public release; distribution unlimited. 


- SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
Monterey, California 93940 


= per." 9 Orr ee Se ae ee 





P13. ABSTRACT 


A model is constructed of a repairable item inventory system where on-hand 
stocks are resupplied both by repair and procurement. The basic properties 
of the model are analyzed. A Poisson failure process is assumed, and steady 
State probabilities are found for operation of the model under <R,Q> and | 
<S-1,S> inventory policies. Further, a compound Poisson failure process under 

an <S-1,S> inventory policy is also analyzed. Measures of supply performance, 
such as expected number of backorders per unit time, are derived for the 


system under the various policies. 





DD hor"..1473 (Pace 1) : 
S/N 0101-807-6811 0 pm Security Clessification = = ae, 








Security Classification 







LINK 8 


KEY WORDS 


Compound Poisson Process 





Inventory Theory 
Poisson Process 
Repairable Item 


Stochastic Process 


————— 


o AE 





7 meer —s 





DO 1473 (Back) 


JN . = 9 DUM MEET GMEEICISoSo oe ee eee ~ 
ec 652) 4) Security Classification Pee ako 




















